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University  of  Texas  at  El  Paso 
Department  of  Geological  Sciences 
Annual  Technical  Report 
Grant  Number  F4962a-J-0438 

Lithospheric  Profiles  in  the  Southwestern  U.S.  Using  Nevada  Test  Site  Sources 
Principal  Investigators:  G.  Randy  Keller,  Diane  I.  Doser,  Kate  C.  Miller 

SUMMARY 

The  objectives  of  this  project  are  to  conduct  a  variety  of  seismic  studies  in  the  southwestern 
U.S.  and  to  use  the  results  to  investigate  lithospheric  structure  and  seismic  wave  propagation  in  the 
area.  In  this  first  year  of  this  project,  good  progress  has  occurred  on  all  fronts.  We  successfully  re¬ 
corded  the  Aiinor  Uncle  explosion  on  ihe  White  Sands  Missile  Range  (WSMR)  in  June  of  1993. 
These  data  have  already  been  played  back ,  edited,  merged,  and  plotted  in  record  section  form.  A 
preliminary  interpretation  is  also  complete.  This  is  a  major  step  in  our  ongoing  effort  to  work  on  a 
lithospheric  transect  extending  between  NTS  and  WSMR. 

We  have  also  made  good  progress  on  the  tomographic  study  of  the  Mojave  desert  area.  Data 
from  the  Cal  Tech  and  U.S.  Geological  Survey  networks  in  the  region  have  been  obtained.  We  have 
implemented  software  for  tomographic  analysis  and  are  close  to  making  initial  runs  on  the  combined 
data  set.  We  gained  valuable  experience  with  tomographic  techniques  by  applying  them  to  Pg  phase 
data  from  the  1989  PACE  experiment.  The  analysis  of  the  large  PACE  data  set  is  nearly  complete. 
Inverse  methods  were  applied  to  the  Pg  phase  from  the  portion  of  the  data  set  which  lies  along  the 
NTS-WSMR  transect  which,  when  combined  with  the  tomographic  results  and  analysis  of  gravity 
anomalies,  provides  an  excellent  picture  on  upper  crustal  structure  in  the  area.  Inversion  of  wide- 
angle  reflections  from  deeper  horizons  in  the  crust  and  upper  mantle  is  almost  complete  and  initial 

results  provide  an  explanation  for  some  of  the  apparently  contradictory  previous  results. 
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The  analysis  of  the  eastern  portion  of  the  NTS-WSMR  profile  was  completed  as  the  subject 
of  an  M.S.  thesis.  Seismic  data  were  recorded  to  offsets  of  1250  km.  These  data  reveal  that  the  deep 
lithospheric  effects  of  the  Rio  Grande  rift  are  more  pronounced  than  predicted  by  regional  models  of 
western  U.  S.  seismic  structure.  In  addition,  the  eastern  extent  of  these  effects  extent  beyond  the 
physiographic  boundaries  of  the  rift  by  about  100  km. 

Our  receiver  function  studies  are  continuing  with  our  analysis  of  the  Lajitas  station  in  west 
Texas  being  almost  complete.  The  regional  analysis  of  surface  wave  dispersion  is  also  nearly  com¬ 
plete. 

The  unexpected  cessation  of  testing  at  NTS  has  caused  some  reorganization  of  our  seismic 
profiling  efforts.  The  conventional  explosion  of  the  Non-Proliferation  Experiment  at  the  NTS,  orig¬ 
inally  planned  for  early  1993,  has  provided  a  very  appropriate  alternative.  This  explosion  was 
delayed  until  September,  1993  because  of  safety  concerns.  This  delay  turned  out  to  be  very  fortu¬ 
itous  because  it  provided  an  opportunity  to  interface  with  a  major  seismic  experiment  targeting  the 
Sierra  Nevada  Mountains.  By  cooperating  with  this  experiment,  we  will  be  able  to  record  a  long 
profile  extending  westward  from  NTS  across  this  mountain  range  and  significantly  add  to  the  overall 
success  of  this  large  cooperative  effort.  A  much  improved  picture  of  earth  structure  around  NTS  will 
result,  and  the  effects  of  the  deep  structure  of  mountain  ranges  on  seismic  waves  is  a  major  concern  in 
efforts  to  monitor  seismic  events  in  many  areas  around  the  world.  This  effort  represents  only  a  minor 
change  in  our  planned  activities  in  that  we  will  be  recording  along  a  profile  extending  westward  from 
NTS  instead  of  southeastward.  We  plan  to  record  along  the  southeastward  profile  in  the  second  year 
of  the  project 

Detailed  descriptions  of  our  main  results  to  date  follow  this  overview  of  our  activities 
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PERSONNEL 

The  Principal  Investigators  on  this  project  are  G.  Randy  Keller,  Diane  I.  Dose,  and  Kate  C. 
Miller.  In  addition,  the  following  students  have  been  involved  to  a  significant  extent: 

Donald  Roberts,  Carlos  Montana,  Donald  Adams,  Bashir  Durrani,  Raed  Aldouri,  Julia  Whitelaw, 
and  Alejandro  Duran.  Donald  Roberts  completed  his  M.  S.  thesis  in  May,  1993  and  Bashir  Durrani 
and  Raed  Aldouri  are  nearing  completion  of  their  Ph.D.  theses.  Mr.  Durrani’s  defense  is  scheduled 
for  September  7. 

Roberts,  Donald  G.,  1993,  A  Long  Offset  Seismic  Profile  Across  the  Rio  Grande  Rift  and  Its  Transi- 
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tion  Zones  in  Central  New  Mexico;  M.  S.  Thesis,  64p. 
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Detailed  Descriptions  of  Specific  Projects 


Tomographic  and  Surface  Wave  Studies  of  the  Mojave  Region 

We  have  compiled  two  important  data  sets  for  studies  of  the  crustal  structure  of  the  Mojave 
Block  and  adjacent  regions,  P-  and  S-wave  arrival  times  for  earthquakes  and  nuclear  tests  within  a 
region  extending  from  114'’W  to  118.5°  W  and  the  IJ.S-Mexico  border  to  38°N  (Figure  1)  recorded 
at  seismograph  stations  operated  by  Caltech,  the  U.S.G.S.  and  the  University  Nevada-Reno  have 
been  collected  for  a  tomographic  study  of  the  region’s  criistal  structure.  We  are  currently  using  avail¬ 
able  geologic  and  geophysical  information  to  construct  a  series  of  block  models  that  will  be  used  as 
starting  models  for  the  tomographic  inversions.  Adiagram  illustrating  the  seismic  ray  coverage  (and 
hence  the  anticipated  resolution  of  velocity  structure  within  the  region)  for  a  selected  group  of  events 
is  shown  in  Figure  2.  Ms.  Julie  Whitelaw  is  conducting  this  study  as  part  of  her  dissertation  research 
under  the  direction  of  Dr.  Doser. 

A  second  data  set  we  have  compiled  is  surface  wave  information  from  teleseisms  located 
along  great  circle  paths  between  the  seismograph  stations  shown  in  Figure  3.  These  data  will  be  used 
to  determine  the  shear  wave  structure  across  the  study  area.  Of  particular  interest  is  the  deep  regional 
shear  velocity  information  obtained  from  surface  waves  with  periods  >  50  sec.  Drs.  Charles  Ammon 
(Univ.  Calif.-Santa  Cruz)  and  Howard  Patton  (Lawrence  Livermore  Lab.)  are  currently  conducting 
a  study  of  surface  waves  from  regional  earthquakes  (periods  <50  sec)  throughout  the  western  United 
States.  We  hope  to  integrate  our  results  form  the  teleseismic  studies  with  their  ongoing  work  to  pro¬ 
duce  a  detailed  model  of  the  shear  velocity  structure  from  the  upper  crust  to  upper  mantle.  Some  of 
the  data  for  our  study  are  available  in  digital  format,  however  seismograms  from  stations  JAS,  ISA 
GSC,  TUC,  and  ALQ  (analog  station  prior  to  ANMO)  must  be  copied  from  microfiche  and  then 
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Figure  1 :  Study  area  for  tomo^aphic  analysis  is  shown  by  bold  line.  P  and  S  arrivals  have  been 
collected  from  seismic  events  (circles)  that  were  recorded  at  the  seismograph  stations  indi¬ 
cated  by  triangles. 
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Figure  2;  Example  of  seismic  ray  path  coverage  from  several  selected  (magnitude  <  5.0)  seismic 
event  within  the  study  area. 
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Figure  3:  Seismograph  stations  (triangles)  to  be  used  in  a  surface  wave  study  to  determine  re¬ 
gional  variations  in  shear  wave  velocity.  Great  circle  paths  between  stations  of  interest  are 
shown  by  lines.  The  dashed  line  between  ANMO  and  TUC  is  a  path  that  has  already  been 
analyzed. 


i  % 


hand  digitized.  Two  undergraduate  studenss,  Ms.  Claudia  Campos  and  Ms.  Anabel  Garcia,  are  work 
ing  under  the  direction  of  Drs.  Doser  and  Keller  on  this  portion  of  the  research. 
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Sierra  Nevada  Foothills  Velocity  Structure  Determined  from  Local 

Earthquake  Network  Data 


Funding  from  this  grant  has  enabled  Dr.  Kate  Miller  to  complete  an  earthquake  study  in  the 
Foothills  of  the  Sierra  Nevada  (Figure  4).  Results  of  this  work  provide  an  important  framework  for 
the  interpretation  of  data  from  the  east-west  transect  of  the  upcoming  Southern  Sierra  Nevada  Conti¬ 
nental  Dynamics  Project  (SSCD)  refraction  experiment.  Her  study  area  lies  immediately  north  of 
the  SSCD  transect.  A  summary  of  these  Sierra  Nevada  Foothills  results  follow.  A  manuscript  on 
these  results  is  currently  in  revision  for  the  Journal  cf  Geophysical  Research  (Miller  and  Mooney, 
1993).  The  Foothills  Metamorphic  Belt  (Figure  4)  is  the  long  band  of  metamorphic  rocks  of  Meso¬ 
zoic  and  Paleozoic  age  that  separates  the  Great  Valley  fore  arc  basin  from  the  Sierra  Nevada  biathlon 
in  northern  and  central  California.  The  belt  is  structurally  complex,  and  has  been  alternatively  inter¬ 
preted  as  a  Jurassic  collisional  zone  or  an  intra-arc  to  forearc  transtensional-transpressional  shear 
system.  Since  the  Foothills  metamorphic  rocks  are  also  the  country  rock  into  which  the  Sierra  Neva¬ 
da  batholith  has  intruded,  late  Mesozoic  plutons  are  also  important  elements  of  the  crustal  structure 
at  the  southern  end  of  the  Foothills  Metamorphic  belt  (Figure  4).  The  deep  crustal  structure  and  com¬ 
position  of  the  Foothills  belt  have  not  been  extensively  study,  but  are  important  to  a  complete  under¬ 
standing  of  the  nature  of  the  deep  crust  in  central  California.  The  presence  of  a  microearthquake 
network  and  a  seismic  reflection  profile  at  the  southern  end  of  the  Foothills  belt  (Figure  4)  permit  the 
first  real  opportunity  to  simultaneously  interpret  the  composition  and  structure  of  the  crust  within 
this  terrane.  Earthquake  hypocenters  occur  at  depths  as  great  as  8  to  40  km  (Figure  5)  which  makes 
the  data  set  particularly  useful  for  obtaining  deepv-crustal  velocity  information  from  seismic  inver¬ 
sion  techniques.  In  this  study,  a  five-layer  velocity  model  is  determined  from  microearthquake  trav¬ 
el-time  data.  The  velocity  model  shows  velocities  of  5.2  to  6.3  km/s  for  the  upper  12  km  of  the  crust, 
and  6.7  to  6.8  km/s  from  12  km  to  an  estimated  Moho  at  32  km  (Figure  6).  The  upper-crustal  veloci- 
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upper  Cenozoic  5©dim©nt5 


Mesozoic  Great  Valley  Sequence 


Mesozoic  Franciscan  Assemblage 


Mesozoic  granitic  rocks 


Mesozoic  metamorphic  rocks  and 
Paleozoic  ophioliDc  rocXs 


Paleozoic  metamorphic  rocks 


I  Figure  4:  Geologic  map  of  study  area.  Chiefrockunitsafter  Jennings  [1977].  Abbrevi- 
1  ations:  MFZ- Melones  Fault  Zone;  BMFZ -  - Bear  Mountains  Fault  Zone;  SAF  — 

San  Andrea.s  Fault;  GIC  -  Guadalupe  Igneous  Complex;  WRP  -  White  Rock  Plu- 
ton;  BLT  -  Bass  Lake  Tonalite.  Solid  triangles  are  earthquake  station  locations. 
Thick  solid  line  is  .seismic  reflection  profile  CC-2.  Dashed  line  in  the  Great  Valley 
is  the  axis  of  the  steepest  gradient  on  the  northeast  side  of  the  Great  Valley  magnetic 
high.  Dot-dash  lin'’  are  refraction  profiles:  B&E  -  Bateman  and  Eaton;  C&M - 
Colburn  and  Mooney;  H&M  -  Holbrook  and  Mooney 
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Figure  5:  A  schematic  cross-section  along  the  CC-2  transect.  The  sub-Foothi!ls  reflec¬ 
tion  pair  are  interpreted  as  shear  zones  formed  during  the  late  Jurassic  Nevadan 
orogeny  that  were  subsequently  truncated  by  one  of  three  mechanisms;  1)  intrusion 
of  batholithic  rocks;  2)  strike-slip  faultingor;  3)  reverse  motion  along  the  near-ver¬ 
tical  shear  zone  that  corresponds  to  the  Bear  Mountains  fault  zone.  Note  that  earth¬ 
quake  epicenters  occur  at  depths  comparable  to  the  shear  zones. 
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Velocity  (km/s) 


Figure  6:  Representative  velocity-dr  h  functions  from  refraction  sur\'eys  and  this  study. 
Velocity  models  discussed  in  the  text  are:  A)  Holbrook  and  Mooney  [1987];  Col¬ 
burn  andMooney  [1986];  B)  Spieth  et  al.  [1981];  C)  This  Study;  D)  Bateman  and 
Eaton  [1967].  Near  10  km  depth  in  the  Holbrook  and  Mooney  model  both  the  op- 
hiolite  and  non-ophiolite  velocitie.s  are  .shown. 
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ties  correspond  to  Foothills  metamorphic  rocks  and  serpentinites  as  well  as  to  diorites  and  granodio- 
rites  of  the  Sierran  batholith,  while  the  lower-cnistal  velocities  are  interpreted  to  represent  inter¬ 
mediate  to  ma;  ^ranulites.  The  velocity  results  differ  significantly  from  other  Foothills  results  and 
from  the  Great  Valley  Province  to  the  west  (Figure  6).  The  low  average  velocities  found  in  the  upper 
12  km  of  the  1-D  velocity  model  are  consistent  with  the  complete  disappearance  of  a  high  velocity 
layer  seen  in  the  Foothills  and  in  the  Great  Valley.  In  its  place,  presumably,  are  variably  serpenti- 
nized  peridotite  with  infolds  of  schists  and  metagraywackes  of  the  Foothills  Belt,  and  granites  and 
granodiorites  of  the  Sierran  bathol  ith.  In  the  lower  crust,  velocities  obtained  in  this  study  are  compa¬ 
rable  to  those  found  beneath  the  Sierra  Nevada  batholith,  but  are  significantly  lower  than  those  seen 
beneath  the  adjacent  Great  Valley.The  majority  of  the  earthquake  hypocenters  as  well  as  a  6.7  kra/s 
layer  in  the  velocity  model  correspond  in  depth  to  thick  zones  of  west-dipping  mid-crustal  reflec¬ 
tions  seen  on  seismic  reflection  profile  CC-2,  that  may  represent  major  shear  zones  formed  during 
the  late  Jurassic  Nevadan  orogeny  or  syn-batholithic  ductile  shear  zones  that  accommodated  crustal 
extension  associated  with  batholith  intrusion.  These  reflections  are  truncated  updip  by  an  inferred 
sub-vertical  contact  that  coincides  with  the  western  edge  of  the  Sierra  Nevada  batholith  and  the 
southward  trace  of  the  Bear  Mountains  fault  zone.  Figure  5  is  a  schematic  cross-section  showing  this 
interpretation  with  earthquake  distribution. 
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Upper  Crustal  Structure  of  Northern  Arizona  Using  Inversion  of  Pg 
Travel-Times  From  the  1989  Pacific  To  Arizona  Crustal  Experiment  (PACE) 

Experiment  Description 

In  order  to  obtain  a  better  understanding  of  the  tectonic  history  and  crustal  structure  of  the  CP, 
a  seismic  study  completely  confined  within  the  CP  was  felt  necessary.  Therefore  during  the  PACE 
experiment  (Figure  7),  data  was  collected  along  an  approximately  230  km  long  crossline  profile 
located  within  the  CP.  This  crossline  profile  extended  in  the  NW-SE  orientation  from  the  north  rim 
of  the  Grand  Canyon  to  northeast  of  Flagstaff,  Arizona  (Figure  7).  This  profile  included  7  small 
in-line  and  3  large  off-end  shots.  These  shots  were  recorded  by  different  portable  instrument  types; 
including  vertical  component  seismic  cassette  recorders  (SCR’s),  vertical  component  seismic  group 
recorders  (SGR’s),  and  vertical  component  PRS-l’s. 

Out  of  470  instruments,  350  were  placed  at  100  m  spacing  centered  about  the  tie  point  of  the 
mainline  (Plateau)  and  crossline  (Grand  Canyon)  profiles.  The  7  small  in-line  shots  were  evenly 
spaced,  about  5  km  apart,  across  this  distance.  The  est  of  the  470  instruments,  total  of  120,  were 
deployed  at  500  m  spacing  to  extend  the  profile  northwest  and  southeast.  This  geometry'  allowed 
recording  of  detailed  wide-angle  reflection  data  from  the  entire  crust  and  possibly  the  upper  mantle. 
Because  of  the  closely  spaced  shots  and  dense  receiver  array,  the  experiment  provided  high  resolu¬ 
tion  data. 

PACE-1989  Data  Description 

The  refraction  data  recorded  during  PACE-1989  along  all  the  profiles  are  of  high  re.solution 
throughout  the  study  area.  This  was  achieved  by  close  spacing  of  both  shots  and  receivers.  Since  the 
development  of  any  seismic  model  depends  on  the  identification  of  various  phases,  a  brief  descrip¬ 
tion  of  the  data  and  various  travel  time  phases  is  given  here. 
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Figure  7:  Map  showing  the  location  of  PACE-1989  mainline  (Plateau)  and  crossline  (Grand 
Canyon)  profiles.  Volcanic  centers  are  shaded. 


16 


Refractions  through  the  highly  fractured,  low  velocity,  sedimentary  rocks  are  observed  up  to 
offsets  of  6.0  km  on  all  the  shots  recorded.  The  apparent  velocities  of  these  arrivals  are  highly  vari¬ 
able  and  range  between  2.0  km/sec  to  more  than  5.0  km/sec.  Thickness  of  this  layer  was  modeled 
only  close  to  the  shot  points  and  interpolated  to  farther  distances. 

The  upper  crustal  Pg  arrivals  which  bottom  in  the  Precambrian  basement  typically  have  ve¬ 
locities  of  6.0  km/sec  and  6.2  km/sec.  The  Pg  refraction  phase  is  very  strong  between  offsets  of  >  6.0 
km  and  100.0  km  after  which  these  are  either  not  observed  or  obscured  by  high  frequency  noise.  The 
Pg  velocity  increases  to  more  than  6.3  km/sec  in  areas  of  mafic  intrusions  into  the  upper  crust. 

A  total  of  1504  Pg  travel  time  picks  were  made  from  the  seismic  sections  up  to  offsets  of  more 
than  100  km.  No  Pg  arrivals  were  picked  from  the  southeastern  most  shot  point  82  because  the  near¬ 
est  receiver  was  more  than  74  km  away  from  the  shot  point  and  the  Pg  phase  could  not  be  clearly 
recognized.  ThePg  travel  time  picks  were  easily  recognizable  on  the  rest  of  the  .seismic  sections  and 
were  of  high  quality.  Tue  Pg  arrivals  w'ere  assigned  an  estimated  average  picking  error  of  0.05  se¬ 
conds. 

Ti'avel-Time  Inversion 

Until  recently,  the  interpretation  of  seismic  refraction/wide  angle  reflection  data  involved 
only  a  tedious  and  painstaking  trial  and  error  ray-trace  forward  modeling  approach.  This  approach 
requires  that  the  model’s  theoretical  response  be  compared  repeatedly  with  the  observed  data  until  a 
velocity  model  is  obtained  that  provides  a  suitable  match  between  the  observed  and  the  calculated 
response.  At  best,  this  requires  countless  person  hours  because  of  many  iterations  required  to  reach  a 
suitable  model.  On  the  other  hand,  seismic  inversion  techniques  are  relatively  fast  and  can  provide 
estimates  of  model  uncertainty,  resolution,  and  non-uniqueness,  but  inversion  techniques  are  not 
always  robust  and  the  parameterization  of  velocity  structure  is  restricted  to  simple  geometries. 
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The  inversion  technique  used  for  this  study  was  developed  by  Zelt  and  Smith  (1992,  Geophs. 
J.  Inti.)  and  was  used  by  them  for  similar  crustal  studies.  This  inversion  technique  is  iterative  and  is 
based  on  a  model  parameterization  and  a  method  of  ray  tracing  suited  to  seismic  refraction/wide 
angle  reflection  data.  The  velocities  and  interfaces  are  defined  as  specified  nodal  points  to  suit  the 
shot-receiver  geometry,  subsurface  ray  coverage,  and  details  of  the  near  surface  structure.  The  sta¬ 
bility  of  inversion  is  improved  by  applying  a  smooth  layer  boundary  simulation  to  a  blocky  mode! 
parameterization.  Rays  are  traced  through  the  velocity  model  using  zero  order  asymptotic  ray  theory 
by  solving  the  ray  tracing  equations  numerically.  The  partial  derivatives  of  travel  time  with  respect 
to  velocity  and  the  depth  of  boundary  nodes  are  calculated  analytically  during  ray  tracing.  Rays  are 
traced  through  each  group  of  arrivals  and  travel  times  and  partial  derivatives  are  interpreted  across 
the  two  closest  ray  end  points.  The  parameter  values  selected  for  adjustment,  both  velocities  and 
boundary  nodes,  are  simultaneously  updated  using  a  damped  least-squares  inversion.  The  inversion 
process  is  stopped  when  either  the  final  model  is  achieved  with  the  desired  trade-off  between  Root- 
Mean-Square  (RMS)  travel-time  residual  and  parameter  residual,  or  rays  are  traced  to  all  observa¬ 
tions.  The  inversion  method  is  particularly  suited  to  multishot  crustal  seismic  refraction/wide  angle 
reflection  data  which  would  otherwise  be  forsvard  modelled. 

Seismic  record  sections  of  4  of  the  10  shot  points  along  the  Grand  Canyon  profile  are  showm 
in  Figures  8,  9, 10,  and  11.  These  data  are  for  shot  points  S3,  85,  76,  and  77,  respectively. 

The  inversion  of  the  PACE-1989  data  was  carried  out  in  four  successive  steps.  First,  the 
upper  crustal  Pg  data  was  inverted  for  velocity  and  upper  crustal  thickness.  Second,  the  mid-crustal 
reflection  data,  PiP  phase,  was  inverted  while  keeping  the  upper  crustal  velocities  and  thickness 
fixed.  The  same  procedure  was  adapted  in  inverting  PmP  phase,  reflections  off  of  the  Moho  and  PtP 
phase  from  within  the  mantle.  This  procedure  was  adapted  to  speed  up  the  inversion  algorithm  and 
obtain  error  estimates  of  the  velocities  and  depths  of  individual  layers.  Since  the  velocities  and 
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Figure  8:  Seismic  record  showing  the  upper  crustal  Pg  arrivals  for  the  crossline  shot  point  83. 

The  vertical  axis  is  reduced  travel-time  with  reduction  velocity  of  6.0  km/sec  while  the 
horizontal  axis  is  shot-to-receiver  offset  in  killometers.  The  traces  have  been  plotted  with  a 
maximum  trace  amplitude  normalization.  The  panels  are  plotted  v/ith  variable  offsets  due  to 
variable  receiver  spacing  in  the  field.  The  upper  crustal  Pg  arrivals  forshot  points  85, 76,  and  77  are 
shown  in  Figures  9,  10,  and  11,  respectively. 


19 


Distance  (km)  *<-  shot  point 


Figure  9:  Seismic  section  for  shot  point  85.  See  Figure  8  for  details. 
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Figure  10;  Seismic  record  for  shot  point  76.  See  Figure  S  for  details. 


21 


NW 


Grand  Canyon 


SE 


^  14 


o 


shot  point  77 


Distance  (km) 


20 

shot  point 


NW 


SE 


Figure  11:  Seismic  record  for  shot  point  77.  See  Figure  8  for  details. 
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depths  of  a  bottom  layer  has  no  effect  on  the  velocity  and  depth  of  the  layer  above  it,  this  procedure 
was  found  most  useful.  A  flat-layer,  constant  velocity  starting  model  was  used  in  each  of  the  steps 
described  above.  The  earth  structure  obtained  in  step  1  was  held  fixed  in  step  2  and  the  structure  in 
steps  1  and  2  was  held  constant  in  step  3  and  so  on.  Since  the  surface  topography  in  the  study  area  was 
added  in  the  model  as  the  top  most  interface,  no  corrections  were  made  to  the  data  for  surface  topog¬ 
raphy  or  near-surface  velocity  variations.  The  elevations  and  depths  in  the  seismic  model  are  with 
respect  to  the  mean  sea  level  with  positive  values  below  sea  level. 

The  cumulative  inversion  result  of  the  upper  crusial  Pg  arrivals  for  the  entire  crossline  is 
shown  in  Figure  12.  For  reasons  of  clarity,  the  model  derived  from  the  inversion  is  plotted  for  indi¬ 
vidual  shots.  Examples  included  are  for  shot  points  83, 85,  and  76,  respectively  (Figures  13, 14,  and 
15). 

In  general,  the  thickness  and  velocities  of  the  final  model  are  well  resolved  with  a  Root  Moan 
Square  (RMS)  error  close  to  the  error  assigned  to  the  travel  time  picLs  (0.069).  The  travel  lime  fit  is 
good  (X'  =  0.211)  and  rays  were  traced  to  almost  all  arrivals. 

Discussion 

The  inversion  result  for  the  upper  crust  can  be  subdivided  into  two  separate  layers  with  velo¬ 
cities  varying  from  2.0  km/sec  to  more  than  5.0  km/sec  in  the  thin,  less  than  2.0  km  thick  sedimenta¬ 
ry  layer  and  about  6.0  kra/s  in  the  Precambrian  basement. 

The  upper  low  velocity,  sedimentary  layer  is  of  variable  thickness  with  thicknesses  varying 
from  less  than  0.5  km  towards  southeast  to  about  1.5  km  towards  northwest.  This  layer  is  observed 
on  all  shots  and  contributes  significantly  to  the  travel  time  delays  observ'ed  in  the  data. 

The  Precambrian  basement  (upper  crust)  could  be  modeled  in  detail  to  depth  of  about  10  kra 
and  contained  significant  velocity  variations.  The  velocity  at  the  upper  boundary  of  this  layer  varies 
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Figure  12:  Cumulative  result  of  the  upper  crustal  velocity  model  across  the  PACE-1989  Grand 
Canyon  profile  derived  from  Pg  arrivals  using  Zell  inversion  algorithm. 

Upper  panel  shows  the  velocity  structure  of  the  upper  10  k.m  of  the  crust  while  the  lower 
panel  shows  the  match  between  observed  (vertical  bars)  and  calculated  (solid  line).  The  data  in 
lower  panel  is  plotted  with  reduction  velocity  of  6.0  km/sec. 
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Figure  13:  Upper  crustal  velocity  model  for  shot  point  83  across  the  PACE-1989  Grand  Canyon 
profile  derived  from  Pg  arrivals  using  Zelt  inversion  algorithm. 
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Figure  14:  Upper  crustal  velocity  model  for  shot  point  85  across  the  PACE-1989  Grand  Canyon 
profile  derived  from  Pg  arrivals  using  2^it  inversion  algorithm. 
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Figure  IS;  Upper  crustal  velocity  model  for  shot  point  76  aaoss  the  PACE-1989  Grand  Canyon 
profile  derived  from  Pg  arrivals  using  Zelt  inversion  algorithm. 
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from  5.5  km/sec  to  about  6.0  km/sec  while  the  velocity  at  the  lower  boundary  varies  from  6.1  kin/sec 
to  about  6.4  km/sec.  The  rays  bottoming  in  this  layer,  the  Pg  arrivals,  are  observed  from  offset  of 
5-10  km  to  offsets  of  greater  then  100  km. 

The  inversion  results  are  in  good  agreement  with  the  geologic  data  in  the  area.  The  velocity 
variations  observed  in  the  Precambrian  basement  can  be  interpreted  in  a  geological  context.  The  low 
velocity  zones  are  interpreted  as  areas  where  the  Late  Proterozoic  Chuar  Group  makes  up  the  upper 
portion  of  the  basement.  This  clastic  unit  outcrops  in  the  Grand  Canyon  and  has  been  encountered  in 
several  wells  in  the  area.  The  high  velocity  zone  under  shot  point  76  could  be  due  to  the  presence  of 
the  Late  Cenozoic  San  Francisco  volcanic  field.  This  area  is  also  seen  as  a  gravity  high  on  Bouguer 
anomaly  maps.  The  velocity  values  obtained  also  correlate  well  with  the  velocities  obtained  by  anal¬ 
ysis  of  sonic  logs  in  wells  in  the  SFVF  and  surrounding  areas. 


Three-Dimensional  Velocity  Structure  of  the  Upper  Crust  Under  the  San 

Francisco  Volcanic  Field,  Arizona,  From  Pg  Arrival  Tomography 

Description  of  the  PACE-1989  Parallel  Seismic  Profile 

The  1989  Pacific  to  Arizona  Crustal  Experiment  (PACE)  included  two  seismic 
refraction/wide  angle  reflection  profiles  crossing  each  other  near  Flagstaff,  Arizona).  Two  shorter 
profiles  (Figure  7)  were  also  recorded  to  provide  some  control  on  basement  velocities.  Profile  3  was 
deployed  by  the  University  of  Texas  at  El  Paso  group  and  consisted  of  an  array  of  22  single 
component  seismic  instruments  in  a  line  parallel  to  the  PACE-1989  mainline.  These  instruments 
recorded  all  the  shots  along  the  PACE  mainline,  and  the  seismic  energy  crisscrossed  through  the 
SFVF.  Tw'o  off-end  shots  were  also  recorded  on  the  receivers  along  the  mainline  and  produced 
seismic  waves  that  also  traveled  through  the  SP/F.  The  ray  coverage  through  the  SFVF  (Figure  1 6) 
provides  an  excellent  opportunity  to  employ  three-dimensional  tomographic  techniques  to  study 
the  seismic  properties  of  the  field.  The  purpose  of  this  study  is  to  investigate  the  possibility  of  an 
upper  crustal  magma  chamber  beneath  San  Francisco  Mountain  and  O’Leary  Peak  and  the  evidence 
for  a  region  of  increased  partial  melt  in  the  upper  crust  beneath  the  volcanic  field  that  would  serve  as 
a  source  for  the  numerous  basaltic  vents. 

The  Seismic  Tomography  Technique 

Tomography  is  defined  as  the  method  of  reconstruction  of  an  object  from  its  projections. 
This  problem  was  first  considered  in  detail  by  Radon  in  1917  in  a  landmark  paper.  Inverse  tech¬ 
niques  based  on  the  Radon  transform  dre  particularly  useful  for  equally  spaced  data  sets,  such  as 
controlled  seismic  refraction  experiments,  since  fast  Fourier  transform  techniques  can  be  imple¬ 
mented.  Radon’s  theorems  can  be  directly  applied  to  uniquely  invert  for  the  image  of  the  object 
through  which  the  specific  rays  are  traveling.  In  seismology,  the  travel  times  can  be  represented  as 
the  line  integral  of  the  slowness  (inverse  of  velocity)  profiles  along  the  specific  ray  paths.  Tomogra- 
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phy  has  found  applications  in  a  wide  variety  of  scientific  fields  such  as  medical  (X-ray,  ultrasonic 
velocity,  or  gamma  ray  intensity),  astrophysics  (radio  astronomy  data),  planetary  science  (occulta- 
tion  data),  and  geophysics  (velocity  inversion).  The  inversion  of  travel  times  into  structure  is  the 
problem  of  imaging  an  object  from  iLs  projections.  In  this  study,  the  time  term  method  of  seismic 
refraction  is  used  to  tomographically  determine  the  slowness  image  and  refractor  topography  of  the 
upper  crust  under  the  San  Francisco  volcanic  field. 

The  time  term  method  of  seismic  refraction  was  first  introduced  in  the  1950’s.  Tho.mas 
Hearn  has  recently  developed  an  expanded  form  of  the  time  term  method  which  could  be  used  to 
determine  delay  times,  lateral  velocity  variations  and  regional  anisotropy,  and  we  have  worked 
closely  with  him  to  implement  his  technique  in  this  study.  The  time  term  method  separates  the 
residual  travel  time  of  a  refracted  ray  into  three  parts  corresponding  to  the  source  delay,  the  station 
delay,  and  the  horizontal  distance  between  source  i  and  station  j.  The  equation  giving  the  travel  time 
of  a  seismic  wave  through  a  structure  consisting  of  two  materials  separated  by  a  plane  interface 
(Figure  17)  may  be  written  in  the  form; 


tij  =  aj  +  bj  +  ~ 


(1) 


where  a;  and  bj  are  ’’time  terms”,  (also  referred  to  as  delay  times  or  time  intercepts)  for  the  sources 
and  receivers,  respectively,  is  the  distance  between  the  shot-point  and  the  seismograph,  tjjisthe 
time  of  propagation  of  a  refracted  wave,  and  v  is  the  velocity  of  uniform  underlying  ’’marker  bed”. 
Equation  (1)  is  linearized  by  using  the  slowness  instead  of  its  reciprocal  velocity,  and  can  be  written 


as: 


tij  =  a;  + 


+  Aij 


(2) 
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Figure  17:  Propagation  of  a  seismic  wave  in  a  two-layer  medium. 
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where  S  is  the  refractor  slowr'ess.  The  station  and  source  delays  could  be  expressed  in  terms  of 


refractor  depth  h,  refractor  slowness  S,  and  the  slowness  depth  profile,  S(z). 
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If  m  recording  stations  observe  waves  from  n  sources,  up  to  n  x  m  equations  of  type  (1)  may 
be  obtained.  The  problem  is  then  to  determine  the  best  value  of  the  velocity  v  and  the  time  terms  for 
all  sources  and  stations.  Finally,  the  time  terms  can  be  interpreted  in  terms  of  the  corresponding 
depths  h,  so  that  the  contours  of  the  interface  can  be  obtained.  The  advantages  of  the  method  are  that 
the  equations  can  be  solved  without  requiring  the  shot  points  and  stations  to  be  laid  out  in  any 
particular  pattern,  the  maximum  amount  of  information  is  extracted  from  the  data,  and  the 
requirements  for  making  simplifying  assumptions  about  the  geological  structure  (assumptions 
which  are,  for  example,  involved  in  describing  the  structure  in  terms  of  plane  layers). 

The  area  under  investigation  is  divided  into  a  grid  of  cells  of  equal  length  and  width.  Each 
ray  is  traced  across  the  grid  and  weighted  sums  of  delays  and  slownesses  are  kept  for  each  cell, 
station,  and  event  involved  with  that  ray.  After  all  of  the  rays  are  used,  the  mean  de’ays  and 
slownesses  are  calculated.  This  process  is  called  back  projection,  since  each  travel  time  is  projected 
back  along  its  ray  path.  Because  data  are  traced  sequentially,  there  is  no  limit  on  the  amount  of  data 
that  can  be  used.  The  computed  delays  depend  on  both  the  crustal  velocity  profile  and  the  refractor 
slowness,  while  the  station  delays  depend  on  both  the  depth  of  the  refractor  and  the  velocity  profile 
above  it.  An  independent  knowledge  of  crustal  velocities  is  desirable  to  separate  these  two  effects. 
However,  the  delay  estimates  can  be  corrected  for  variations  in  refractor  velocity  by  assuming  a 
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constant  offset  distance,  the  horizontal  distance  that  the  ray  travels  from  the  refractor  to  the  station. 
This  offset  distance  can  be  calculated  by  expanding  equation  3  by  Taylor’s  series  which  gives; 


Si  (sq)  =  Si  (Sq  +  ^  ^^0  Si  (Sq) 

s  h 

3;  (Sq)  3;  (^0  +  Ss)  +  i==6s 

"  s- 

(sq)  =  (sq  +  6s)  +  F  6s 
where 

F  =  =  13.0  bn 

The  parameter  F  is  dependent  upon  the  unknown  velocity  profile,  therefore,  it  is  approximated  by  a 
constant  value  of  13  km.  This  value  assumes  an  average  upper  crustal  velocity  of  5.5  km/sec,  a 
refractor  velocity  of  6.0  km/sec,  and  a  refractor  depth  of  5  km.  The  velocity  and  depth  estimates  of 
the  upper  crustal  layer  and  the  refractor  depth  used  here  as  input  were  obtained  from 
two-dimensional  inversion  of  the  data  along  the  crossline.  All  station  statics  are  corrected  to  this 
constant  offset  value.  The  corrected  delays  obtained  are  interpreted  as  either  refractor  depth  or 
crustal  velocity  variations  (slownesses).  In  the  case  of  a  constant  velocity  crust  of  Vq  =  5.5  km/sec, 
and  refractor  velocity  of  vj  =  6.0  km/sec,  an  a^  =  0.1  sec  of  relative  delay  corresponds  to: 


h 


=  1.4  km 
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of  refractor  topography.  In  the  case  of  a  flat  reflector,  the  same  delay  will  correspond  to  a  0.5  km/sec 
upper  crustal  velocity  change. 

The  estimated  slownesses  will  be  interpreted  to  delineate  areas  of  positive  and  negative 
delays  associated  with  an  increase  or  decrease  in  the  crustal  \eiocities. 

Tomographic  Imaging  of  the  SFVF  Using  Pg  Arrivals 

Due  to  the  short  shot-receiver  offset  along  the  parallel  line,  the  majority  of  the  arrivals  pass¬ 
ing  through  the  SFVF  were  identified  as  upper  crustal  Pg  arrivals.  In  this  study,  data  processing  was 
started  by  picking  Pg  arrivals  from  the  travel-time  data.  This  amounted  to  more  than  4,000  arrivals. 
Pg  phase  picks  between  13  km  and  160  km  were  windowed  and  plotted  on  a  time-distance  plot  and  a 
straight  line  was  fitted  to  the  data.  The  outliers  of  more  than  1  sec  from  the  straight  line  fit  were  re¬ 
moved.  The  data  in  Figure  18  are  plotted  with  the  reduction  velocity  of  6.0  km/sec,  which  is  the  in¬ 
verse  slope  of  the  line  fitted  to  the  data.  The  tomographic  process  is  iterative.  In  the  first  step,  the 
delays,  slowness  image,  and  the  scaling  factors  (one  to  multiply  all  the  event  delays,  one  to  multiply 
all  the  station  delays,  and  one  to  multiply  the  slowness  image)  were  calculated.  After  this,  new  resid¬ 
uals  were  calculated.  These  residuals  were  then  back  projected,  and  the  resulting  slowness  perturba¬ 
tions  were  added  to  the  previous  estimates.  This  process  was  repeated  until  the  change  in  the  Root 
Mean  Square  (RMS)  error  was  very  small.  All  results  shown  here  are  for  fifteen  iterations.  The  slow¬ 
ness  results  are  shown  in  Figure  19,  and  the  associated  delays  are  shown  in  Figure  20.  The  mean 
intercept  and  the  inverse  siopte  of  the  line  fitted  to  the  data  in  Figure  18  is  0.64  sec  and  6.0  km/sec, 
respectively.  The  mean  velocity  perturbations  in  Figure  19  are  ±  0.25  km/sec.  This  velocity  range 
agrees  well  with  the  velocities  obtained  by  two-dimensional  inverse  modelling  of  the  data  set  along 
the  crossline  profile.  The  mean  station  delays  are  very  small  (0.003  sec).  Negative  delays,  denoting 
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Figure  18:  Pg  travel-time  pick  are  plotted  as  a  function  of  source  to  receiver  offset. 
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eariy  arrivals,  are  observed  towards  the  edge  of  the  parallel  line  closer  to  the  southeastern  boundary 
of  the  SFVR 


The  quality  of  the  solution  obtained  is  determined  by  inverting  artificial  travel  times  and 
then  inverting  a  synthetic  model.  This  allows  ns  to  empirically  investigate  the  resolution  and  vari¬ 
ances  of  the  slowness  image  obtained  for  the  actual  data.  By  applying  the  tomographic  algorithm  to 
these  aritificial  travel  times,  the  variance  of  the  velocity  values  determined  can  be  estimated.  Instead 
of  actual  arrival  times,  travel  times  that  have  a  Gaussian  distribution  witli  a  standard  deviation  of 
0.055  sec  assigned  to  the  highest  quality  arrivals,  were  used.  This  is  a  large  amount  of  error.  TIte 
actual  ray  paths  of  the  dataset  were  used  in  the  reconstruction.  The  result  of  inverting  these  synthetic 
data  indicate  that  Figure  r>  is  a  reliable  depiction  of  the  upper  crustal  velocity  variations  in  the  area. 

Discussion 

The  most  significant  feature  of  the  tomography  result  in  Figure  19  is  the  zone  of  high  com- 
pressional  velocity  at  upper  crustal  depths  (less  than  10  km)  below  San  Francisco  Mountain.  The 
higher  upper  cnistal  Pg  velocities  of  6.3  to  6.9  km/sec  correlate  well  with  intermediate  to  basaltic 
rocks  of  the  San  Francisco  volcanic  field.  The  higher  upper  crustal  velocities  suggest  the  intrusion  of 
magmas  of  basaltic  composition  into  the  upper  crust  below  San  Francisco  Mountain  stratovolcano  as 
the  volcano  formed.  If  most  of  the  dikes  and  sills  that  have  intruded  the  upper  crust  in  the  SFV'F  have 
cooled  below  their  solidus  temperatures,  the  average  velocity  of  the  intruded  region  would  be  higher 
than  the  highly  fractured  and  moresilicicbasement  rocks  surrounding  it.  The  high  upper  crustal  vel- 
icity  of  6.7  km/sec  to  6.9  km/sec  under  the  San  Francisco  Mountain  compare  well  with  experimen¬ 
tally  derived  P-wave  velocit'es  for  gabbros  at  upper  crustal  pressures.  Gabbro,  the  intrusive  equiva¬ 
lent  of  basalt,  has  P-wave  velocities  of  6.83  km/scc  near  the  surface  to  about  6.99  km/sec  at  depths  of 
30  km.  The  high  velocity  regions  also  correlate  well  with  positive  gravity  anomalies. 
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The  time  term  analysis  technique  provides  reasonable  velocity  structure  of  the  San  Francisco 
Volcanic  field.  High  velocity  zones  are  detected  under  San  Francisco  Mountain,  and  O’lxary  Peak 
extending  within  4  km  of  surface  to  about  7  km  in  depth.  The  P-wave  velocities  within  these  zones 
are  from  5%  to  14%  higher  than  the  surrounding  basement  rocks.  These  high  velocities  could  he 
caused  by  the  intrusion  of  numerous  mafic  dikes  and  sills  into  upper  crust,  cooling  to  subsol  idus  tem¬ 
peratures  while  the  volcanic  centers  grew.  Gravity  modeling  of  the  San  Franci.sco  Volcanic  field  in¬ 
dicates  that  the  high  velocity  regions  are  coincident  with  the  high  density  bodies  at  upper  crustal 
depths  under  many  of  the  volcanic  centers.  It  could,  therefore,  be  said  with  confidence  that  the  high 
velocities  found  here  are  evident  in  various  data  sets  that  have  been  interpreted  independently  of 
each  other. 
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Long  Offset  Profile  Across  the  Rio  Grande  Rift  and  its  Transition  Zones  in 

Central  New  Mexico 

A  number  of  geophysical  surveys  have  been  conducted  over  the  last  20  years  to  examine  the 
intra-rift  structure  and  morphology  of  the  Rio  Grande  Rift  (RGR).  The  purpose  of  this  research  is  to 
expand  on  these  previous  studies  by  using  long  offset  seismic  data  to  examine  the  transition  zone 
between  the  southern  RGR  and  the  neighboring  Great  Plains  (GP)  province  to  the  east  The  profile 
of  recordings  begins  southeast  of  Socorro,  New  Mexico  and  extends  eastward  to  near  the  Texas- 
New  Mexico  border  (Figure  21).  The  surface  expression  of  rifting  within  this  region  is  postulated  to 
extend  as  far  east  as  the  western  extreme  of  the  Sacramento  and  Pedemal  Uplifts  but  no  seismic  in¬ 
vestigations  have  examined  the  transition  at  depth. 

Tlie  profile  was  chosen  so  it  would  have  a  near  radial  orientation  to  the  Nevada  Test  Site 
(NTS).  The  large  nuclear  explosions  created  at  NTS  provide  excellent  seism  ic  sources  and  the  resul¬ 
tant  long  offsets  (  >  1000  km)  make  it  possible  to  examine  the  structure  of  the  transition  region  to 
depths  below  the  low  velocity  zone  in  the  upper  mantle.  The  investigation  targeted  these  depths  be¬ 
cause  it  is  presumed  that  the  RGR-GP  transition  extends  into  the  deep  lithosphere  and 
asthenosphere. 

However,  these  long  offsets  also  present  a  problem.  Since  the  ultimate  objective  of  the  sur¬ 
vey  was  to  analyze  the  rift/craton  transition  zone,  a  modeling  technique  that  allowed  for  two-dimen¬ 
sional  velocity  structures  was  necessary.  The  method  chosen  for  this  survey  was  a  2-D  ray  tracing 
algorithm  based  on  asymptotic  ray  theory  (ART).  In  addition,  since  all  of  the  station  offsets  in  this 
work  were  over  900  km,  some  type  of  compensation  had  to  be  made  for  the  curvature  of  the  earth. 
Conventional  earth  “flattening”  transforms  involve  various  mathematical  techniques  for  modifying 
a  homogeneous,  flat  earth  model  so  that  the  results  obtained  from  the  specific  modeling  algorithm 
approximate  those  expected  from  a  homogeneous,  spherical  earth.  The  two-dimensional  nature  of 
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the  models  used  in  ART  make  the  application  of  such  a  transformation  difficult.  To  alleviate  this 
problem,  a  new  technique  of  compensating  for  the  curvature  was  developed  for  this  study. 

We  implemented  a  geometric  solution  to  the  curvature  problem:  the  earth  was  approximated 
by  a  series  of  pie-shaped  wedges.  Within  each  of  these  wedges,  the  earth  was  assumed  to  be  flat,  but 
the  cumulative  error  in  this  assumption  is  small  if  the  pie  wedges  are  kept  small,  e.g.  about  50  kra.  A 
true  Cartesian  coordinate  system  is  set  up  with  the  X-axis  tangent  to  the  earth’s  surface  at  the  center 
ofthe  user’s  model  and  with  Y  increasing  towards  the  center  of  the  earth.  Each  of  the  interface.s  in  the 
user’s  “flat-earth”  model  space  is  divided  into  a  series  of  line  segments,  and  then  each  of  the  end¬ 
points  of  the  line  segments  is  transformed  from  the  user’s  coordinate  system  into  the  true  Cartesian 
system  using  the  mathematical  relationship  described  in  Figure  22.  The  result  is  a  series  of  curv-ed 
surfaces  drawn  in  a  truly  Cartesian  coordinate  system. 

The  modeling  algorithm  then  proceeds,  using  the  pC,  Y)  coordinate  system,  to  determine  the 
appropriate  ray  paths  and  travel  times.  Each  ray  is  stopped  at  the  boundary  corresponding  to  the 
earth’s  surface  (as  opposed  to  Y=0)  and  the  final  (X,  Y)  coordinates  are  inverse  transformed  to  deter¬ 
mine  the  range  at  which  the  ray  exited  the  model.  The  accuracy  of  this  technique  was  tested  by 
comparing  the  its  output  with  that  predicted  by  conventional  earth  flattening  algorithms  using  a  ho¬ 
mogeneous,  layered  earth  model.  Careful  observation  of  these  data  show  that  this  approach  accu¬ 
rately  predicts  the  arrival  times  to  distances  of  greater  than  1000  km. 

Analysis  of  Seismic  Sections 

The  seismic  sections  gathered  during  this  work  are  displayed  in  Figure  23.  The  first  arrivals 
across  both  sections  are  interpreted  as  critical  refractions  from  the  Moho  discontinuity.  The  se¬ 
cond  arrival  is  Pi:  an  arrival  refracted  below  the  upper  mantle  low  velocity  zone.  Synthetic  sections 
generated  using  a  standard  earth  model  for  the  region  suggest  that  several  phases  with  aixival  times  of 
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Figure  22:  Graphic  representation  of  earth  curvature  algorithm. 

The  user  provides  (R,  D)  of  point  P  and  the  algorithm  returns  (Xp,  Yp). 

P  =  subsurface  position  of  some  point  within  the  user’s  model 
E  =  point  on  surface  directly  above  point  P 
S  =  origin  for  user’s  range/depth  coordinate  system 
R  =  arc  length  from  S  to  E  (user  input  —  range) 

D  =  radial  distance  from  E  to  P  (user  input  —  depth) 

O  =  origin  of  Cartesian  coordinate  system 

M  =  range  at  which  the  true  Cartesian  coordinate  system  is  tangential  to  the 
curved  earth 
A  =  center  of  the  earth 
H  =  radius  of  the  earth 
C  =  arc  distance  from  S  to  M 
Xm  =  distance  from  O  to  M  (=  C) 

9=  (R-C)/H 
Xp=  Xn.  +  (H-D)*sin0 
Yp=  H-(H-D)*cos9 


44 


Time  -  Di.stance/8.0  (km) 


1000  1050  1100  1150  1200  1250 


I  Distance  (km) 


Figure  23:  Seismic  sections  from  NTS2  and  NTS4  with  phase  arrivals  labeled. 

The  top  section  is  from  shotpoints  NTSl  and  NTS2  and  the  bottom  section  is  from 
NTS4.  The  inset  shows  a  closeup  of  the  arrivals  from  NTS4. 
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about  10  sec  are  converging  within  the  aperture  of  this  work  (1000  km  to  1200  km).  This  is  due 
partly  to  the  nature  of  the  synthetic  model  (several  thin  layers  at  depth)  but  the  observed  sections 
indicate  similar  multiple  arrivals.  Further  analysis  of  these  tightly  clustered  phases  will  have  to  wait 
for  additional  traces  to  be  added  to  this  survey  in  order  to  increase  the  aperture  and  make  phase  cor¬ 
relations  easier.  It  should  be  noted  however  that  in  the  observed  sections,  this  packet  of  energy  and 
the  Pn  arrivals  arrive  significantly  later  (~2  sec)  than  predicted  by  the  ID  models.  In  addition,  the  PI 
arrival  does  not  appear  to  be  as  continuous  across  both  sections  as  the  P^  arrival.  A  significant  delay 
in  the  PI  arrival  is  apparent  between  offsets  of  R  =  1050  and  R  =  1150. 

Initial  Model 

The  initial  model  generated  for  this  survey  was  derived  from  several  sources.  Crustal  models 
were  chosen  from  surveys  within  each  of  the  separate  geologic  provinces  that  exist  between  NTS  and 
the  eastern  extremes  of  the  seismic  array  (Basin  and  Range,  Colorado  Plateau,  Rio  Grande  Rift,  and 
Great  Plains)  and  an  upper  mantle  model  was  drawn  from  regional  studies  in  the  western  and  central 
U.S.  The  combined  model  consisted  of  three  crustal  layers  and  several  mantle  layers  extending  to  a 
depth  of  about  200  km. 

Table  1  summarizes  the  various  iterations  used  to  model  the  Pn  arrival.  The  best  fit  to  the  data 
is  obtained  with  an  “abrupt,”  50-lan  wide  transition  zone  between  rift  and  cratonic  crust  extending 
from  R  =  1150  km  to  R=  1200  km.  The  crust  thickens  from  32  km  to  51  km  yielding  an  apparent  dip 
along  the  Moho  of  approximately  21°.  This  centers  the  transition  directly  under  the  Pecos  River. 
Bringing  the  transition  closer  to  the  rift  (farther  west)  caused  the  modeled  phases  to  arrive  as  much  as 
1  second  later  than  the  observed  phases.  Similarly,  moving  the  transition  to  the  east  caused  the  mod¬ 
eled  rays  to  arrive  up  to  1  sec  too  early.  Analysis  of  a  gravity  profile  along  the  trace  of  this  seismic 
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line  shows  a  large  gravity  gradient  in  eastern  New  Mexico  which  correlates  with  the  general  location 
of  the  crustal  transition  derived  from  this  seismic  study. 


Table  1:  Description  of  model  iterations  on  the  Moho  transition  within  the  study  area.  The  highlighted 

box  represents  the  best  fit. 


Rl 

Moba . . 1  Rt 

. 

Rl 

R2 

Description 

Results 

1050 

1200 

Broad  transition  beginning  under 
the  Sacramentos  and  extending  to 
eastern  banks  of  Pecos  River. 

Arrivals  up  to  R=1150  km  are 
very  good.  East  of  R=1150,  arriv¬ 
als  are  late,  especially  between 
R=1175  and  R=1225  (>0.5  sec 
late,  up  to  1.2  sec  late). 

1050 

1100 

Extreme  case  with  abrupt  transi¬ 
tion  to  cratonic  crust  under  Sacra¬ 
mento  Mountains. 

Western  arrivals  OK.  Model  trav¬ 
el  times  begin  to  arrive  too  late  by 
R=1120.  By  R=1200,  predicted 
arrival  i.<;  ~1  sec  late. 

1100 

1150 

Abrupt  transition  beginning  at 
eastern  end  of  Capitan  Intrusion 
and  ending  25  km  west  of  Pecos 
River. 

All  arrivals  east  of  R=1190  are 
about  V2  sec  late. 

1150 

1200 

Abrupt  transition  centered  on  the 
Pecos  River  and  extending 
approx.  25  km  to  either  side. 

Good  fit  across  all  Pn  arrivals. 
Arrivals  between  1100  and  1150 
are  very  slightly  early  (0.1  -  0.15 
sec). 

1250 

1260 

Extreme  case  where  no  seismic 
energy  reaches  the  transition. 

Most  arrivals  from  R=1100  km 
eastward  are  too  early  (0.1  -  1.0 
sec). 

Fewer  constraints  were  available  for  the  deeper  parts  of  the  model.  As  a  first  guess,  the 
Burdick  and  Helmberger  (1978)  T7  model  was  used  across  this  entire  model  beneath  65  km.  The 
resultant  PI  arrivals  were  about  1  sec  early  across  most  of  the  section  and  over  2  sec  early  in  the 
region  ofR=1100km.  However,  by  lowering  the  velocity  in  each  of  the  layers  under  65  km  depth  by 
only  0.06  km/sec,  it  was  possible  to  get  a  reasonable  fit  across  most  of  the  section.  The  only  mis- 
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match  is  in  the  region  of  the  apparent  delay  in  the  observed  PI  arrival  time  around  R  =  1100  km.  On 
closer  inspection,  a  very  weak  arrival  can  be  detected  across  this  entire  region  which  is  co  ainuous 
with  the  larger  arrival  across  the  rest  of  the  section.  Since  the  data  in  each  trace  is  normali2ed  to  a 
fixed  amplitude  height,  a  stronger  arrival  just  behind  the  first  PI  arrival  dominates  the  signal  in  that 
time  window.  Therefore,  the  phase  match  for  PI  is  consistent  across  the  entire  section. 

The  exact  source  of  the  large  secondary  arrival  in  the  PI  window  is  impossible  to  determine 
with  this  small  data  set.  However,  one  possibility  is  energy  that  initially  entered  the  bottom  of  the 
Colorado  Plateau  upper  mantle  but  was  then  refracted  by  the  large  velocity  variation  in  the  upper 
mantle  between  the  Plateau  and  the  rift  eventually  reached  the  surface  within  the  rift.  It  was  po.ssible 
to  get  some  arrivals  to  do  this  through  the  model  and  they  did  have  the  proper  arrival  times,  but  the 
solution  was  very  unconstrained  and  non-unique. 

Figure  24a  is  a  drawing  of  the  final  model  derived  from  this  work.  The  figure  was  drawn  as  a 
“flat  earth,”  i.e.  without  the  earth  curvature  approximations  discussed  above.  Figure  24b  shows  all 
the  rays  which  were  used  to  construct  the  final  model.  Since  the  data  had  to  be  divided  into  2  separate 
sections,  the  modeling  was  also  divided.  The  solid  rays  are  those  that  originated  at  shotpoint 
NTS1/NTS2  and  the  dashed  rays  are  those  that  originated  at  NTS4.  Figure  25  shows  the  NTS2  and 
NTS4  seismic  sections  plotted  side-by-side  and  overlaid  by  the  modeled  arrival  branches  and  Fig¬ 
ure  26  uses  velocity  vs.  depth  profiles  to  show  the  differences  between  this  model  within  the  Rio 
Grande  Rift  and  Burdick  and  Helmberger’s  (1978)  model  T7. 

CONCLUSIONS 

A  number  of  conclusions  can  be  drawn  from  the  unreversed  data  presented  here.  First,  the 
location  of  the  crustal  transition  between  the  Rio  Grande  Rift  and  the  Great  Plains  is  about  50  km  or 
more  farther  east  than  any  surface  features  attributed  to  the  rift  (Figure  27).  This  indicates  that  the 
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Figure  25;  NTSl/NTS2and  NTS4  seismic  section  overlaid  by  modeled  arrival  branches. 

The  station  at  1149  km  was  occupied  for  both  shots.  Data  to  its  left  is  from 
NTS1/NTS2  and  the  rest  is  from  NTS4. 
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Figure  2^?:  Velocity /Depth  profiles  within  the  Rio  Grande  Rift  from  the  final  model  pres¬ 
ented  here  and  the  Burdick  and  Helmberger  (1978)  model  T7. 

The  solid  line  represents  the  T7  model  and  the  dashed  line  is  the  model  from  this 
study.  The  profile  was  taken  at  R  =  1000  km. 
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slight  off-line  location  of  the  shotpoints. 


crustal  thinning  and  associated  mantle  upwarp  within  the  rift  have  undercut  the  westernmost  Great 
Plains  Province.  It  further  confirms  that  the  Tularosa  basin  is  well  within,  and  a  product  of,  the  Rio 
Grande  Rift. 

By  comparing  the  synthetic  sections  generated  from  the  ID  model  of  Burdick  and 
Helmberger  with  the  observed  sections  it  is  clear  that  the  recorded  energy  in  the  study  area  traveled 
through  significantly  slower  material  than  that  predicted  by  the  T7  model.  Furthermore,  the  rift  crust 
in  the  final  model  of  this  work  is  actually  thinner  and  faster  than  that  in  the  Burdick  and  Helmberger 
model  (Figure  26).  Therefore,  the  compressional  velocities  within  the  upper  mantle  and/or  astheno- 
spheric  structure  of  the  modeled  area  must  be  slower  than  the  respective  T7  structures.  Unfortu- 
neately,  the  resolution  of  this  survey  is  not  sufficient  to  identify  the  exact  location  of  the  low  velocity 
region.  The  anomaly  can  be  accounted  for  with  a  0.06  km/s  adjustment  to  the  T7  structures  across  the 
entire  modeled  area.  This  results  in  a  less  than  1%  drop  in  velocity  across  the  lithosphere/astheno- 
sphere  transition.  However,  the  anomalously  high  surface  heat  flow  and  voluminous  volcanics  with¬ 
in  the  rift  suggests  that  the  low  velocity  anomaly  may  be  isolated  under  the  rift  and  would  therefore 
represent  more  than  just  a  0.06  km/s  velocity  drop.  Teleseismic  investigations  of  the  rift  have  indi¬ 
cated  a  depth  to  the  low  velocity  zone  of  about  70  km  under  the  north  central  part  of  the  rift  with  an 
8%  velocity  reduction.  Burdick  and  Helmberger  (1978)  place  the  interface  at  a  depth  of  65  km  for  an 
average  of  the  western  U.S.  with  a  velocity  drop  of  3.7%.  Even  using  the  upper  mantle  velocities 
from  Burdick  and  Helmberger  (1978),  a  velocity  drop  of  8%  in  the  upper  asthenosphere  beneath  the 
rift  would  result  in  velocities  (~7.4  km/s)  that  are  intuitively  too  low  for  that  depth  and  would  not 
match  the  refraction  data  presented  here.  The  investigators  who  conducted  the  teleseismic  study 
acknowledge  that  the  8%  velocity  reduction  is  extremely  large  but  indicate  that  such  a  value  is  neces¬ 
sary  to  produce  the  large  teleseismic  delays  they  observed.  The  discrepancy  between  the  teleseismic 
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results  and  refraction  results  may  be  explained  by  the  different  locations  of  the  two  profiles  or  the 
different  paths  that  the  seismic  energy  traveled  before  reaching  the  respective  stations. 
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Prelimmary  Results  from  the  Refraction  Profile  Extending  Westward  from  the 
Minor  Uncle  Conventional  Explosion  at  White  Sands  Missile  Range 

As  part  of  our  ongoing  investigation  of  crustal  structure  between  the  White  Sands  Missile 
Range  (WSMR)  and  the  Nevada  Test  Site  (NTS),  UTEP  deployed  a  seismic  array  in  June.  1993,  to 
record  the  Minor  Uncle  conventional  surface  explosion  at  WSMR.  Three  recording  crews  deployed 
a  total  of  68  stations  at  1  mile  intervals  along  US  Highway  60  beginning  west  of  Socorro,  NM,  and 
ending  several  miles  west  of  Datil,  NM  (Figure  28).  Each  station  consisted  of  a  1  Hz,  vertical  com¬ 
ponent,  Mark  Products  L-4  seismometer  and  a  USGS  Seismic  Cassette  Recorder  (SCR).  These  SCR 
units  are  automatic  gain  ranging,  analog  tape  recorders  built  by  the  USGS.  We  used  these  older  units 
because  of  the  scheduling  difficulties  associated  with  the  newer  Seismic  Group  Recorders  (SGRs)  or 
the  PASSCAL  Reftek  instruments.  Our  previous  experiences  with  the  SCRs  have  convinced  us  of 
their  high  reliability  and  excellent  recording  quality.  The  only  limitation  of  the  SCRs,  a  narrow 
recording  window,  was  handled  by  maintaining  good  communications  between  the  recording  crew's 
and  the  personnel  at  WSMR.  In  addition  to  the  68  SCRs,  five  digital  Reftek  recorders  with  three- 
component  seismometers  were  deployed  at  widely  spaced  intervals  along  the  line.  These  modem 
instruments  were  deployed  primarily  for  us  to  acquire  more  experience  with  them  and  to  compare 
their  recording  characteristics  with  those  of  the  SCRs. 

The  Minor  Uncle  experiment  went  off  as  scheduled  and  66  SCRs  and  all  five  Refteks  re¬ 
corded  a  strong  signal.  The  resultant  record  section  (Figure  29)  'S  extremely  sharp  and  shows  a  num¬ 
ber  of  prominent  arrivals.  The  PmP  phase  marked  on  the  section  represents  the  wide-angle  reflec¬ 
tions  off  of  the  crust/raantle  transition  and  is  a  strong  indicator  of  overall  crustal  structure.  The 
approximately  7  sec  arrival  time  of  this  phase  at  about  100  kra  compares  well  with  much  lower  reso¬ 
lution  surveys  conducted  just  to  the  south  of  this  profile.  At  the  farthest  offsets,  it  appears  that  the 
PmP  signal  is  merging  with  Pq,  the  refracted  arrival  from  the  upper  mantle.  The  most  remarkable 


55 


Figure  28;  Basemap  showing  the  location  of  the  Minor  Uncle  conventional  surface  explosion  and 
the  seismic  stations  deployed  by  UTEP. 
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arrival  on  this  section  is  the  late  arrival  that  crosses  the  western  half  of  the  section  at  about  13  se¬ 


conds.  This  appears  to  be  a  reflection  from  within  the  upper  mantle.  Work  is  in  progress  to  fully 
interpret  these  data  and  to  evaluate  their  implications  for  the  NTS-WSMR  lithospheric  transect. 
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